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ABSTRACT 

An  integrally  stiffened  graphite/epoxy  composite 
rotorcraft  structure  is  evaluated  via  computational 
simulation.  A computer  code  that  scales  up  con- 
stituent micromechanics  level  material  properties 
to  the  structure  level  and  accounts  for  all  possible 
failure  modes  is  used  for  the  simulation  of  compos- 
ite degradation  under  loading.  Damage  initiation, 
growth,  accumulation,  and  propagation  to  fracture 
are  included  in  the  simulation.  Design  implications 
with  regard  to  defect  and  damage  tolerance  of  inte- 
grally stiffened  composite  structures  are  examined. 
A procedure  is  outlined  regarding  the  use  of  this 
type  of  information  for  setting  quality  acceptance 
criteria,  design  allowables,  damage  tolerance,  and 
retirement-for-cause  criteria. 

KEYWORDS:  Composites,  Composite  Shells, 
Composite  Structures,  Computational  Simulation, 
Damage,  Degradation,  Durability,  Fracture,  Inte- 
grated Stiffening,  Laminates,  Structural  Degrada- 
tion. 

INTRODUCTION 

Laminated  composite  structures  are  used  in  many 
aerospace  applications  such  as  rotorcraft  compo- 
nents, advanced  aircraft  fuselage,  rocket  motor 
cases,  pressure  vessels,  containment  structures,  and 
other  components  with  various  shapes  and  sizes. 
In  these  applications  composite  structures  are  re- 
quired to  withstand  significant  in-plane  loads.  Ad- 
ditionally, composite  structures  are  required  to 
possess  sufficient  bending  stiffness  to  resist  buck- 
ling. Discussion  in  the  current  paper  is  focussed 
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on  integrally  stiffened  composite  structures  sub- 
jected to  in-plane  loads  and  internal/external  pres- 
sures. Damage  initiation,  growth,  accumulation, 
and  propagation  to  fracture  is  simulated  for  com- 
posite panels  and  cylindrical  shells  with  and  with- 
out integrated  stiffening  layups.  The  influence  of 
integrated  stiffeners  is  examined  with  regard  to  out 
of  plane  stiffness  contribution  as  well  as  damage 
progression  and  structural  durability  assessment 
under  applied  loading.  Changes  in  the  damage  ini- 
tiation load  and  the  structural  fracture  load  are 
quantified  due  to  the  presence  of  integrated  stiff- 
eners. 

Integrally  stiffened  composite  structures  are  man- 
ufactured by  adding  additional  plies  or  braids  that 
are  designed  to  improve  both  the  out-of-plane  stiff- 
ness and  buckling  resistance  of  composite  structural 
components.  Integrated  stiffening  is  typically  ob- 
tained within  the  manufacturing  process  of  a com- 
posite material.  The  layup  of  a laminated  compos- 
ite or  the  braid  structure  of  a braided  composite 
may  be  spatially  concentrated  to  obtain  a stiffen- 
ing of  the  composite  in  the  desired  orientations.  In- 
tegrally stiffened  composites  would  display  logical 
patterns  of  ply/braid  concentrations.  For  example, 
a rectangular  composite  panel  may  be  stiffened  by 
adding  bands  of  layers  with  fibers  along  the  diago- 
nals of  the  panel.  Another  example  is  the  superpo- 
sition of  lattice  type  periodic  angled  stiffeners  over  a 
smooth  composite  plate  layup  during  the  manufac- 
turing process.  Integrated  stiffening  of  composites 
may  be  obtained  by  hand  layup  of  composite  tape, 
filament  winding,  or  by  automated  braiding  pro- 
cesses. The  distinguishing  feature  of  an  integrally 
stiffened  composite  is  that  the  stiffening  system  has 
been  manufactured  simultaneously  as  part  of  the 
stiffened  composite  structure. 
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Critical  components  of  a structure  are  required  to 
remain  safe  and  be  able  to  function  under  load- 
ing after  experiencing  some  damage.  The  cause 
of  damage  may  be  an  accident,  defect,  or  unex- 
pected overloading.  Damage  tolerance  of  a struc- 
ture is  quantified  by  the  residual  strength,  that  is 
the  additional  load  carrying  ability  after  damage. 
Composite  structures  are  well  suited  for  design  with 
emphasis  on  damage  tolerance  as  continuous  fiber 
composites  have  the  ability  to  arrest  cracks  and  pre- 
vent self-similar  crack  propagation.  It  is  difficult  to 
design  and  certify  a composite  structure  because  of 
the  complexities  in  predicting  the  overall  congruity 
and  performance  of  fiber  composites  under  various 
loading  and  hygrothermal  conditions. 

Design  considerations  with  regard  to  the  durabil- 
ity of  fiber  composite  structures  require  an  a priori 
evaluation  of  damage  initiation  and  fracture  propa- 
gation mechanisms  under  expected  loading  and  ser- 
vice environments.  Concerns  for  safety  and  surviv- 
ability of  critical  components  require  a quantifica- 
tion of  the  structural  fracture  resistance  under  load- 
ing. A significant  design  parameter  with  regard  to 
composite  damage  tolerance  is  the  laminate  config- 
uration. In  general,  quasiisotropic  laminates  yield 
better  damage  tolerance.  However,  in  many  cases 
a quasiisotropic  laminate  may  not  be  the  most  ef- 
ficient with  regard  to  structural  strength  and  per- 
formance when  there  is  no  damage.  For  a ratio- 
nal design  process  it  is  necessary  to  quantify  the 
structural  damage  tolerance  for  a candidate  design. 
The  ability  of  designing  composites  with  numerous 
possible  fiber  orientation  patterns,  choices  of  con- 
stituent material  combinations,  ply  drops  and  hy- 
bridizations, render  a large  number  of  possible  de- 
sign parameters  that  may  be  varied  for  an  optimal 
design.  Damage  initiation  and  progression  charac- 
teristics are  much  more  complex  for  integrally  stiff- 
ened fiber  composites  compared  to  homogenous  or 
orthotropic  materials.  The  structural  fracture  pro- 
cess of  a fiber  composite  depends  on  many  parame- 
ters such  as  laminate  configuration,  fiber  volume 
ratio,  constituent  stiffness/s trength/hygrothermal 
parameters,  stiffening  system,  and  the  fabrication 
process.  Recent  developments  in  computational 
simulation  technology  have  made  it  possible  to  eval- 
uate the  details  of  progressive  damage  and  fracture 
in  composite  structures.  Computational  simulation 
enables  assessment  of  the  damage  initiation  and 
propagation  loads.  A damage  energy  release  rate 
is  evaluated  globally  during  simulation  by  comput- 
ing the  work  done  per  unit  damage  created.  The 
damage  energy  release  rate  is  used  to  quantify  the 
structural  damage  tolerance  at  different  stages  of 
degradation.  The  influence  of  local  defects  or  flaws 


and  effects  of  the  fabrication  process  in  terms  of 
residual  stresses  are  taken  into  account. 

An  important  feature  of  computational  simulation 
is  the  assessment  of  damage  stability  or  damage  tol- 
erance of  a structure  under  loading.  At  any  stage 
of  damage  progression,  if  there  is  a high  level  of 
structural  resistance  to  damage  progression  under 
the  service  loading,  the  structure  is  stable  with  re- 
gard to  fracture.  The  corresponding  state  of  struc- 
tural damage  is  referred  to  as  stable  damage.  On 
the  other  hand,  if  damage  progression  does  not 
encounter  significant  structural  resistance,  it  cor- 
responds to  an  unstable  damage  state.  Unstable 
damage  progression  is  characterized  by  very  large 
increases  in  the  amount  of  damage  due  to  small 
increases  in  loading.  Whereas  during  stable  dam- 
age progression  the  amount  of  increase  in  damage 
is  consistent  with  the  increase  in  loading. 

Internal  damage  in  composites  is  often  initiated  as 
cracking  due  to  normal  stresses  transverse  to  fiber 
orientation.  At  the  presence  of  stress  concentra- 
tions or  defects,  initial  damage  may  also  include 
fiber  fracture.  Further  degradation  is  in  the  form  of 
additional  fiber  fractures  that  usually  lead  to  struc- 
tural fracture.  Because  of  the  numerous  possibili- 
ties with  material  combinations,  composite  geom- 
etry, fiber  orientations,  and  loading  conditions,  it 
is  essential  to  have  an  effective  computational  ca- 
pability to  predict  the  behavior  of  composite  struc- 
tures for  any  loading,  geometry,  composite  mate- 
rial combinations,  and  boundary  conditions.  The 
predictions  of  damage  initiation,  growth,  accumu- 
lation, and  propagation  to  fracture  are  important 
in  evaluating  the  load  carrying  capacity  and  reli- 
ability of  composite  structures.  Quantification  of 
the  structural  fracture  resistance  is  also  required  to 
evaluate  the  durability /life  of  composite  structures. 

Laminated  composite  design  practice  has  been 
based  on  extensive  testing  with  attempts  to  apply 
formal  fracture  mechanics  concepts  to  interpret  test 
results.  In  certain  cases  interpretation  of  laminated 
composite  test  data  via  fracture  mechanics  has  been 
satisfactory.  However,  in  most  cases  fracture  me- 
chanics methods  have  significantly  mispredicted  the 
strength  of  fiber  composites.  Reconciliation  of  test 
results  with  fracture  mechanics  has  required  signif- 
icant modifications  of  effective  fracture  toughness 
and  specific,  laminate  configuration  dependent,  ef- 
fective stress  concentration  field  parameters.  Addi- 
tionally, required  adjustments  of  fracture  mechan- 
ics parameters  have  had  to  be  reassessed  with  every 
change  in  constituent  and  laminate  characteristics. 
The  complete  evaluation  of  laminated  composite 
fracture  requires  an  assesment  of  ply  and  subply 
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level  damage/fracture  processes. 

Continuous  fiber  composites  in  general  have  the 
ability  to  arrest  cracks  and  prevent  self-similar 
crack  propagation.  For  most  fiber  reinforcement 
configurations,  cracks  and  other  stress  concentra- 
tors do  not  have  as  important  an  influence  in  com- 
posites as  they  do  for  homogeneous  materials.  An- 
other important  aspect  is  the  multiplicity  of  design 
options  for  composites.  The  ability  of  designing 
composites  with  numerous  possible  fiber  orienta- 
tion patterns,  choices  of  constituent  material  com- 
binations, ply  drops,  hybridizations,  and  integrated 
stiffening  options  render  a large  number  of  possible 
design  parameters  that  may  be  varied  for  an  opti- 
mal design. 

The  present  approach  by-passes  traditional  frac- 
ture mechanics  to  provide  an  alternative  evaluation 
method,  conveying  to  the  design  engineer  a detailed 
description  of  damage  initiation,  growth,  accumula- 
tion, and  propagation  that  would  take  place  in  the 
process  of  ultimate  fracture  of  an  integrally  stiff- 
ened fiber  composite  structure.  Results  show  in 
detail  the  damage  progression  sequence  and  struc- 
tural fracture  resistance  during  different  degrada- 
tion stages.  This  paper  demonstrates  that  com- 
putational simulation,  with  the  use  of  established 
material  modeling  and  finite  element  modules,  ad- 
equately tracks  the  damage  growth  and  subsequent 
propagation  to  fracture  for  an  integrally  stiffened 
fiber  composite  structure. 

METHODOLOGY 

Computational  simulation  is  implemented  via  the 
integration  of  three  modules:  (1)  composite  me- 
chanics, (2)  finite  element  analysis,  and  (3)  dam- 
age progression  tracking.  The  overall  evaluation  of 
composite  structural  durability  is  carried  out  in  the 
damage  progression  module  (Minnetyan  et  al  1990) 
that  keeps  track  of  composite  degradation  for  the 
entire  structure.  The  damage  progression  module 
relies  on  the  composite  mechanics  code  (Murthy 
and  Chamis  1986)  for  composite  micromechanics, 
macromechanics  and  laminate  analysis,  and  calls  a 
finite  element  analysis  module  that  uses  anisotropic 
thick  shell  elements  to  model  laminated  composites 
(Nakazawa  et  al  1987). 

A computational  simulation  cycle  begins  with  the 
definition  of  constituent  properties  from  a materials 
databank.  Composite  ply  properties  are  computed 
by  the  composite  mechanics  module.  The  compos- 
ite mechanics  module  also  computes  through-the- 
thickness  structural  properties  of  each  laminate. 
The  finite  element  analysis  module  accepts  the  com- 


posite properties  that  are  computed  by  the  compos- 
ite mechanics  module  at  each  node  and  performs 
the  analysis  for  a load  increment.  After  an  incre- 
mental finite  element  analysis,  the  computed  gener- 
alized nodal  force  resultants  and  deformations  are 
supplied  to  the  composite  mechanics  module  that 
evaluates  the  nature  and  amount  of  local  damage,  if 
any,  in  the  plies  of  the  composite  laminate.  Individ- 
ual ply  failure  modes  are  determined  using  failure 
criteria  associated  with  the  negative  and  positive 
limits  of  the  six  ply-stress  components  ( am , <7/22, 
0733,  0/12)  0/i3,  0-/23),  a modified  distortion  energy 
(MDE)  combined  stress  failure  criterion.  The  MDE 
failure  criterion  is  expressed  as: 

*«w5!Sf“-(gs)i  (i) 

where  a and  /?  indicate  tensile  or  compressive 
stress,  5/n  a is  the  local  longitudinal  strength  in 
tension  or  compression,  5/22<*  is  the  transverse 
strengthen  tension  or  compression,  and 

X , — (1  + 4^/12  - vnz)Et 22  + (1  - ^/23)E/n 

^ [E/11E/22  (2  + Vm  + 3 ) (2  + IS£21  "H^m)]1/2 

(2) 

The  type  of  failure  is  assessed  by  comparison  of 
the  magnitudes  of  the  squared  terms  in  Equation 
(1).  Depending  on  the  dominant  term  in  the  MDE 
failure  criterion,  fiber  failure  or  matrix  failure  is 
assigned. 

The  generalized  stress-strain  relationships  for  each 
node  are  revised  according  to  the  composite  dam- 
age evaluated  by  the  composite  mechanics  module 
after  each  finite  element  analysis.  The  model  is 
automatically  updated  with  a new  finite  element 
mesh  and  properties,  and  the  structure  is  reana- 
lyzed for  further  deformation  and  damage.  If  ply 
failure  criteria  indicate  new  or  additional  damage 
during  a load  increment,  the  damage  tracking  mod- 
ule degrades  the  composite  properties  affected  by 
the  damage  and  reanalizes  the  structure  under  the 
same  load.  When  there  is  no  indication  of  fur- 
ther damage  under  a load,  the  structure  is  consid- 
ered to  be  in  equilibrium.  Subsequently,  another 
load  increment  is  applied  leading  to  possible  dam- 
age growth,  accumulation,  or  propagation.  In  the 
computational  simulation  cases  presented  in  this 
paper,  analysis  is  stopped  when  commencement  of 
the  damage  propagation  phase  is  indicated  by  lam- 
inate fracture.  Laminate  fracture  is  predicted  when 
major  principal  failure  criteria  are  met  for  all  plies 
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at  a node.  After  laminate  fracturing,  the  compos- 
ite structure  is  anticipated  to  enter  a final  damage 
propagation  stage  that  leads  to  ultimate  structural 
fracture  or  collapse. 

During  progressive  damage  tracking  the  follow- 
ing terminology  is  utilized  to  describe  the  various 
stages  of  degradation  in  the  composite  structure: 
(1)  damage  initiation  refers  to  the  start  of  damage 
induced  by  loading  that  the  composite  structure  is 
designed  to  carry;  (2)  damage  growth  is  the  pro- 
gression of  damage  from  the  location  of  damage 
initiation  to  other  regions;  (3)  damage  accumula- 
tion is  the  increase  in  the  amount  of  damage  in  the 
damaged  regions  with  additional  damage  modes  be- 
coming active;  (4)  damage  propagation  is  the  rapid 
progression  of  damage  to  other  regions  of  the  struc- 
ture; (5)  structural  fracture  is  the  ultimate  disinte- 
gration of  the  specimen. 

INTEGRALLY  STIFFENED 
PANEL 

A graphite/epoxy  laminated  composite  plate  with 
integrated  ±45°  intermittent  lattice  stiffeners  was 
investigated  with  damage  and  fracture  propagation 
due  to  tension,  compression,  and  in-plane  shear 
loads.  The  response  of  the  integrally  stiffened  com- 
posite panel  was  compared  with  that  of  an  unstiff- 
ened  skin  plate.  The  unstiffened  plate  was  given 
additional  skin  thickness  such  that  the  material  vol- 
ume was  the  same  as  the  material  volume  of  the 
integrally  stiffened  plate.  The  additional  plies  of 
the  unstiffened  plate  were  given  fiber  orientations 
of  ±45°,  same  as  the  fiber  orientations  of  the  in- 
tegrated stiffeners.  Both  the  unstiffened  and  inte- 
grally stiffened  panels  were  made  of  AS-4  graphite 
fibers  in  a high  modulus  high  strength  epoxy  ma- 
trix (AS-4/HMHS).  Ply  layup  of  the  unstiffened 
panel  was  [±45/0/90 / ±45]g . The  stiffened  panel 
had  a skin  layup  of  [0/90/±45]g.  The  stiffener  plies 
were  added  to  the  top  of  the  skin  as  [+45]$  or  as 
[-45]$.  Ply  thickness  was  0.127  mm  (0.005  in).  The 
fiber  volume  ratio  was  Vj=0.60  and  the  void  vol- 
ume ratio  vas  Vt,=0.01.  The  cure  temperature  was 
Tcu=177°C  (350°F).  The  fiber  and  matrix  prop- 
erties were  obtained  from  a databank  of  compos- 
ite constituent  material  properties  resident  in  the 
composite  mechanics  module  (Murthy  and  Chamis 
1986).  The  fiber  and  matrix  properties  correspond- 
ing to  this  case  are  given  in  Tables  I and  II,  respec- 
tively. The  HMHS  matrix  properties  are  represen- 
tative of  the  3501-6  resin.  These  in-situ  properties 
are  similar  to  those  identified  by  Sobel  et  al  (1993) 
with  experimental  correlation. 


The  panels  had  identical  planar  geometry  with  a 
width  of  305  mm  (12.0  in.)  and  length  of  457  mm 
(18  in.).  Each  finite  element  model  contained  117 
nodes  and  96  uniformly  sized  square  elements.  Fig- 
ure 2 shows  the  finite  element  model  with  diagonal 
lines  along  the  ±45°  stiffener  bands.  Numbers  at 
nodes  indicate  the  laminate  type  number  at  that 
node  for  the  integrally  stiffened  panel.  Laminate 
type  1 had  a layup  of  [0/90/±45]g,  representing 
the  skin.  Laminate  type  2 had  the  skin  layup  plus 
[+45]$  representing  the  +45  stiffened  nodes.  Lam- 
inate type  3 had  the  skin  layup  plus  [-45]$  repre- 
senting the  -45  stiffened  nodes.  Laminate  type  4 
had  the  skin  layup  plus  [+45$/-45$]  representing  in- 
tersection nodes  for  +45  and  -45  stiffeners.  Panels 
were  assumed  to  be  simply  supported  for  structural 
response;  i.e.  the  left  end  nodes  of  the  panel  were 
restrained  against  all  displacement  components  and 
the  right  end  nodes  were  restrained  against  dis- 
placement normal  to  the  plane  of  the  panel.  Ad- 
ditionally, the  right  end  nodes  were  constrained  to 
have  uniform  displacement  in  the  plane  of  the  panel 
via  duplicate  node  specifications. 

Structural  response  characteristics  of  the  integrally 
stiffened  panel  were  evaluated  in  terms  of  the  buck- 
ling load  and  the  bending  stiffness.  Analysis  of  the 
integrally  stiffened  panel  under  uniaxial  compres- 
sion indicated  a buckling  load  of  1,114  N (250  lbs). 
The  buckling  load  of  the  unstiffened  panel  was  only 
213  N (48  lbs).  Therefore  the  buckling  resistance  of 
the  integrally  stiffened  panel  was  5.2  times  that  of 
the  unstiffened  panel  with  the  same  amount  of  com- 
posite material.  Similarly,  the  bending  rigidity  was 
significantly  improved  due  to  integrated  stiffeners. 
Figure  3 shows  a comparison  of  the  midspan  de- 
flections of  unstiffened  and  integrally  stiffened  pan- 
els due  to  uniform  bending  moment  applied  to  the 
simply  supported  right  end  of  the  panel.  Figure 
3 shows  that  the  integrally  stiffened  panel  was  7.2 
times  stiffer  than  the  unstiffened  panel  in  bending. 

Next,  the  in-plane  progressive  damage  and  fracture 
responses  were  compared  for  the  unstiffened  and 
stiffened  panels.  Figure  4 shows  damage  progres- 
sions of  unstiffened  and  integrally  stiffened  panels 
subjected  to  tension.  Damage  initiation  for  the  in- 
tegrally stiffened  panel  occurred  in  the  90°  skin 
plies  due  to  transverse  tensile  fractures.  Damage 
initiation  for  the  unstiffened  panels  was  also  in  the 
90°  plies  due  to  transverse  tensile  fractures.  In  both 
cases  the  MDE  combined  stress  failure  criterion  was 
activated.  However,  the  damage  initiation  load  for 
the  integrally  stiffened  panel  was  approximately  one 
third  of  the  damage  initiation  load  for  the  flat  panel 
with  the  same  material  volume.  Figure  5 shows  the 
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displacements  in  tension.  The  uniaxial  stiffnesses 
of  the  unstiffened  and  integrally  stiffened  panels  are 
approximately  the  same  prior  to  damage  initiation. 
However,  stiffness  of  the  integrally  stiffened  panel 
degrades  quickly  due  to  damage  initiation  and  pro- 
gression. Figure  6 shows  damage  progressions  of 
the  unstiffened  and  integrally  stiffened  panels  un- 
der uniaxial  compression.  Damage  initiation  for  the 
integrally  stiffened  panel  occurred  at  the  edge  of  the 
panel  at  midspan  where  +45  and  -45  stiffeners  con- 
verged (laminate  type  4).  The  damage  initiation 
modes  included  the  transverse  tensile  and  longitu- 
dinal compressive  failures  of  the  90°  skin  plies,  in- 
plane shear  failures  of  the  ±45°  skin  plies,  as  well  as 
the  in-plane  shear  failures  of  the  +45  stiffener  plies 
near  the  skin.  Damage  initiation  for  the  unstiffened 
panels  under  compression  was  in  the  0°  plies  due 
to  longitudinal  compressive  fractures.  The  dam- 
age propagation  load  of  the  unstiffened  panel  was 
more  than  five  times  that  of  the  integrally  stiffened 
panel.  The  initial  stiffnesses  of  the  unstiffened  and 
integrally  stiffened  panels  in  compression  were  the 
same  as  observable  from  Figure  7.  However,  the 
stiffness  of  the  integrally  stiffened  panel  degraded 
at  a much  lower  loading  compared  to  degradation 
load  of  the  flat  panel.  Figure  8 shows  the  damage 
progressions  of  unstiffened  and  integrally  stiffened 
panels  under  in-plane  simple  shear  loading.  Dam- 
age initiation  for  the  integrally  stiffened  panel  was 
due  to  longitudinal  compressive  failure  of  a +45 
skin  ply  at  a type  3 stiffened  node.  On  the  other 
hand,  damage  initiation  for  the  flat  panel  under 
shear  was  due  to  transverse  tensile  failures  of  the 
90°  plies.  As  it  was  in  tension  and  compression,  also 
in  shear  the  integrally  stiffened  panel  degraded  at 
a lower  load  compared  to  the  flat  panel.  Figure  9 
shows  the  load-displacement  relationships  in  shear. 
Computational  simulation  results  depicted  in  Fig- 
ure 9 indicate  that  in  shear,  the  initial  stiffness  of 
the  flat  panel  was  slightly  higher  than  that  of  the 
integrally  stiffened  panel. 

INTEGRALLY  STIFFENED 
COMPOSITE  SHELL 

A composite  cylindrical  shell  with  and  without  in- 
tegrated ±45°  intermittent  lattice  stiffeners  was  in- 
vestigated with  damage  and  fracture  propagation 
due  to  axial  tension,  as  well  as  internal  and  exter- 
nal pressure  loads.  The  response  of  the  integrally 
stiffened  shell  was  compared  with  that  of  an  un- 
stiffened shell.  The  unstiffened  shell  was  given  ad- 
ditional skin  thickness  such  that  the  material  vol- 
ume was  the  same  as  the  material  volume  of  the 
integrally  stiffened  shell.  The  additional  plies  of 


the  unstiffened  shell  were  given  fiber  orientations 
of  ±45°,  same  as  the  fiber  orientations  of  the  in- 
tegrated stiffeners.  Both  the  unstiffened  and  in- 
tegrally stiffened  shells  were  made  of  AS-4/HMHS 
composite.  Ply  layup  of  the  unstiffened  shell  was 
[±45/0/90/±45]5 . The  stiffened  shell  had  a skin 
layup  of  [0/90/±45]g.  The  stiffener  plies  were 
added  to  the  top  of  the  skin  as  [+45]g  or  as  [- 
45]g.  Ply  thickness  was  0.127  mm  (0.005  in).  The 
fiber  volume  ratio  was  V/=0.60  and  the  void  vol- 
ume ratio  vas  Vv=0.01.  The  cure  temperature  was 
Tcu=177°C  (350°F). 

Each  cylindrical  shell  finite  element  model  con- 
tained 400  nodes  and  384  uniformly  sized  square 
elements.  Figure  10  shows  the  finite  element  model 
for  a cylindrical  shell.  The  shell  was  simulated  sub- 
ject to  increasing  levels  of  internal  and  external 
pressurizations  as  well  as  axial  loading.  To  rep- 
resent the  axial  stresses  produced  in  the  closed  end 
pressure  vessel,  boundary  nodes  at  one  end  of  the 
cylinder  were  subjected  to  force  resultants  in  the 
axial  direction,  whereas  the  boundary  FEM  nodes 
were  restrained  in  the  axial  direction  at  the  oppo- 
site end.  The  uniformly  distributed  axial  tension 
was  such  that  the  generalized  axial  stresses  in  the 
shell  wall  were  half  those  developed  in  the  hoop  di- 
rection. 

Figure  11  shows  a comparison  of  damage  progres- 
sions for  the  integrally  stiffened  and  unstiffened 
composite  shells  under  axial  tension.  Figure  11  in- 
dicates that  ultimate  strength  of  the  integrally  stiff- 
ened cylindrical  shell  is  approximately  60  percent 
of  the  ultimate  strength  of  shell  without  integral 
stiffeners.  Additionally,  damage  initiation  for  the 
integrally  stiffened  shell  begins  at  only  10  percent 
of  the  ultimate  load.  On  the  other  hand,  damage 
initiation  for  the  unstiffened  shell  corresponds  to  a 
loading  level  that  is  approximately  97  percent  of  its 
ultimate  load.  Figure  12  shows  the  exhausted  cu- 
mulative damage  energy  based  on  depleted  energies 
of  the  local  failure  modes.  The  exhausted  damage 
energy  represents  a similar  pattern  of  damage  pro- 
gression as  the  percent  damage  volume  depicted  in 
Figure  12.  Figure  13  shows  the  damage  energy  re- 
lease rates  (DERR)  based  on  the  incremental  work 
done  by  applied  loading  during  damage  progression 
of  the  integrally  stiffened  cylindrical  shell.  Peak 
values  in  the  DERR  levels  indicate  significant  dam- 
age events.  The  first  peak  in  DERR  correspond- 
ing to  damage  initiation  for  the  integrally  stiffened 
shell  occurred  due  to  transverse  tensile  fractures  in 
the  90°  skin  plies.  The  steep  increase  in  DERR 
at  approximately  11  kN  (2.5  k)  axial  tension  cor- 
responds to  a partial  separation  between  skin  and 
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stiffener  plies. 

Figure  14  shows  comparison  of  damage  progressions 
for  integrally  stiffened  and  unstiffened  cylindrical 
shells  subjected  to  external  pressures.  The  ulti- 
mate pressure  for  the  integrally  stiffened  shell  is 
approximately  75  percent  of  the  ultimate  pressure 
for  the  unstiffened  shell.  Damage  initiation  for  the 
integrally  stiffened  shell  occurs  at  23  percent  of  its 
ultimate  pressure.  On  the  other  hand  damage  ini- 
tiation for  the  unstiffened  shell  occurs  at  approxi- 
mately 94  percent  of  its  ultimate  pressure.  Also,  the 
volume  of  damage  in  the  integrally  stiffened  shell  at 
ultimate  pressure  is  much  greater  than  the  volume 
of  damage  in  the  unstiffened  shell  at  ultimate  pres- 
sure. 

SUMMARY  OF  RESULTS 

This  section  summarizes  some  of  the  insights 
gained  by  the  present  investigation  as  applicable 
to  graphite/epoxy  laminated  composite  structures 
with  integrated  stiffeners,  as  follows: 

1.  Integrated  stiffeners  may  be  used  to  increase 
the  out-of-plane  stiffness  properties  signifi- 
cantly. 

2.  The  improvement  of  out  of  plane  stiffness  is 
paid  for  by  the  degradation  of  in-plane  damage 
tolerance  and  durability  characteristics. 

3.  Damage  initiation  occurred  much  sooner  for 
the  integrally  stiffened  composite  structures 
subjected  to  in-plane  tension,  compression,  or 
shear. 

4.  Damage  initiation  occurred  in  the  skin  plies  for 
all  cases  considered. 

5.  For  the  integrally  stiffened  composites,  dam- 
age initiation  was  within  the  skin  at  stiffened 
nodes. 


GENERALIZATION  OF 
PROCEDURE 

The  present  computational  simulation  method  is 
suitable  for  the  design  and  continued  in-service 
evaluation  of  composite  structures.  Compos- 
ite structures  with  different  constituents  and  ply 
layups  can  be  evaluated  under  any  loading. 

Structural  health  monitoring  is  based  on  damage 
tolerance  requirements  defined  via  the  computa- 
tional simulation  method.  A damage  tolerance  pa- 
rameter is  described  as  the  state  of  damage  after  the 
application  of  a given  load  level,  normalized  with 


respect  to  the  damage  state  corresponding  to  ulti- 
mate fracture.  Identification  of  damage  progression 
mechanisms  and  the  sequence  of  progressive  frac- 
ture modes  conveys  useful  information  to  evaluate 
structural  safety.  Computational  simulation  results 
can  be  formulated  into  health  monitoring  criteria, 
increasing  the  reliability  of  composite  structures. 
The  simulated  failure  modes  and  the  type  of  failure 
provide  the  necessary  quantitative  and  qualitative 
information  to  design  an  effective  health  monitor- 
ing system.  Computed  local  damage  energy  release 
rates  are  correlated  with  the  magnitudes  of  acous- 
tic emission  signals  and  other  damage  monitoring 
means  such  as  piezoelectric  stress  sensors  and  strain 
gages  that  are  an  integral  part  of  a composite  struc- 
ture. Fiber  optics  data  networks  embedded  in  the 
composite  structure  would  transmit  the  detected 
local  damage  information  to  an  expert  system  that 
provides  feedback  and  reduces  power  to  delay  fail- 
ure. 

The  basic  procedure  is  to  simulate  a computational 
model  of  the  composite  structure  subjected  to  the 
expected  loading  environments.  Various  fabrication 
defects  and  accidental  damage  may  be  represented 
at  the  ply  and  constituent  levels,  as  well  as  at  the 
laminate  level.  Computational  simulation  may  be 
used  to  address  various  design  and  health  monitor- 
ing questions  as  follows: 

1.  Evaluation  of  damage  tolerance:  Computa- 

tional simulation  will  generate  the  damage  that 
would  be  caused  due  to  overloading  by  the  type 
of  load  the  structure  is  designed  to  carry.  On 
the  other  hand,  a fabrication  defect  or  acci- 
dental damage  produced  by  inadvertent  load- 
ing that  is  not  an  expected  service  load  can 
be  included  in  the  initial  computational  model. 
Once  the  composite  damage  is  defined,  damage 
tolerance  can  be  evaluated  by  monitoring  dam- 
age growth  and  progression  from  the  damaged 
state  to  ultimate  fracture.  Significant  param- 
eters that  quantify  damage  stability  and  frac- 
ture progression  characteristics  are  the  rate  of 
damage  increase  with  incremental  loading,  and 
the  changes  in  the  structural  response  charac- 
teristics with  loading.  Identification  of  damage 
initiation/progression  mechanisms  and  the  se- 
quence of  progressive  fracture  modes  convey 
serviceable  information  to  help  with  critical 
decisions  in  the  structural  design  and  health 
monitoring  process.  Determination  of  design 
allowables  based  on  damage  tolerance  require- 
ments is  an  inherent  use  of  the  computational 
simulation  results.  Simulation  of  progressive 
fracture  from  defects  allows  setting  of  qual- 
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ity  acceptance  criteria  for  composite  structures 
as  appropriate  for  each  functional  requirement. 
Detailed  information  on  specific  damage  toler- 
ance characteristics  help  establish  criteria  for 
the  retirement  of  a composite  structure  from 
service  for  due  cause. 

2.  Determination  of  sensitive  parameters  affect- 
ing structural  fracture:  Computational  simu- 
lation indicates  the  damage  initiation,  growth, 
and  progression  modes  in  terms  of  a damage 
index  that  is  printed  out  for  the  degraded  plies 
at  each  damaged  node.  In  turn,  the  dam- 
age index  points  out  the  fundamental  physi- 
cal parameters  that  characterize  the  compos- 
ite degradation.  For  instance,  if  the  dam- 
age index  shows  ply  transverse  tensile  fail- 
ure, the  fundamental  physical  parameters  are 
matrix  tensile  strength,  fiber  volume  ratio, 
matrix  modulus,  and  fiber  transverse  modu- 
lus; of  which  the  most  significant  parameter 
is  the  matrix  tensile  strength  (Murthy  and 
Chamis,  1986).  In  addition  to  the  significant 
parameters  pointed  out  by  the  ply  damage  in- 
dex, sensitivity  to  hygrothermal  parameters 
may  be  obtained  by  simulating  the  composite 
structure  at  different  temperatures  and  mois- 
ture contents.  Similarly,  sensitivity  to  resid- 
ual stresses  may  be  assessed  by  simulating  the 
composite  structure  fabricated  at  different  cure 
temperatures.  Identification  of  the  important 
parameters  that  significantly  affect  structural 
performance  for  each  design  case  allows  opti- 
mization of  the  composite  for  best  structural 
performance.  Sensitive  parameters  may  be 
constituent  strength,  stiffness,  laminate  config- 
uration, fabrication  process,  and  environmen- 
tal factors. 

3.  Interpretation  of  experimental  results  for  de- 
sign decisions:  Computational  simulation  al- 
lows interactive  experimental-numerical  as- 
sessment of  composite  structural  performance. 
Simulation  can  be  used  prior  to  testing  to  iden- 
tify locations  and  modes  of  composite  damage 
that  need  be  monitored  by  proper  instrumen- 
tation and  inspection  of  the  composite  struc- 
ture. Interpretation  of  experimental  data  can 
be  significantly  facilitated  by  detailed  informa- 
tion from  computational  simulation.  Subscale 
experimental  results  may  be  extended  to  full 
prototype  structures  without  concern  for  scale 
effects  since  computational  simulation  does  not 
presume  any  global  parameters  but  is  based  on 
constituent  level  damage  tracking. 


CONCLUSIONS 

On  the  basis  of  the  results  obtained  from  the  in- 
vestigated flat  composite  plate,  integrally  stiffened 
panel,  integrally  stiffened  and  unstiffened  compos- 
ite cylindrical  shell  examples  and  from  the  general 
perspective  of  the  available  computational  simula- 
tion method,  the  following  conclusions  are  drawn: 

L Computational  simulation  can  be  used  to  track 
the  details  of  damage  initiation,  growth,  and 
subsequent  propagation  to  fracture  for  unstiff- 
ened and  integrally  stiffened  composite  struc- 
tures. 

2.  For  the  considered  integral  stiffening  system, 
out-of-plane  structural  response  characteristics 
such  as  the  buckling  load  and  the  bending  stiff- 
ness are  significantly  improved. 

3.  In-plane  load  carrying  capability  of  a compos- 
ite panel  is  reduced  due  to  damage  initiation 
and  progression  processes  caused  by  the  pres- 
ence of  integrated  stiffeners. 

4.  Computational  simulation,  with  the  use  of  es- 
tablished composite  mechanics  and  finite  ele- 
ment modules,  can  be  used  to  predict  the  in- 
fluence of  composite  geometry  as  well  as  load- 
ing and  material  properties  on  the  durability 
of  composite  structures. 

5.  The  demonstrated  procedure  is  flexible  and  ap- 
plicable to  all  types  of  constituent  materials, 
structural  geometry,  and  loading.  Hybrid  com- 
posites and  homogeneous  materials,  as  well  as 
laminated,  stitched,  woven,  and  braided  com- 
posites can  be  simulated. 

6.  Computational  simulation  by  CODSTRAN 
represents  a new  global  approach  to  progres- 
sive damage  and  fracture  assessment  for  any 
structure. 
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TABLE  I:  AS-4  Fiber  Properties: 


Number  of  fibers  per  end  = 10000 

Fiber  diameter  = 0.00762  mm  (0.300E-3  in) 

Fiber  Density  = 4.04E-7  Kg/m3  (0.063  lb/in3) 
Longitudinal  normal  modulus  = 200  GPa  (29.0E-F6 
psi) 

Transverse  normal  modulus  = 13.7  GPa  (1.99E+6 
psi) 

Poisson’s  ratio  (^12)  = 0.20 

Poisson’s  ratio  (^23)  = 0.25 

Shear  modulus  (G12)  = 13.8  GPa  (2.00E+6  psi) 

Shear  modulus  (G23)  = 6.90  GPa  (1.00E+6  psi) 

Longitudinal  thermal  expansion  coefficient  = -1.0E- 

6/°C  (-0.55E-6  /°F) 

Transverse  thermal  expansion  coefficient  ~ 1.0E- 
5/°C  (0.56E-5  /°F) 

Longitudinal  heat  conductivity  = 301 

kJ-m/hr/m2/°C  (4.03  BTU-in/hr/in2/°F) 
Transverse  heat 

conductivity  — 30.1  kJ-m/hr/m2/°C  (.403  BTU- 
in/hr/in2/°F) 

Heat  capacity  = 0.712  kJ/kg/°C  (0.17  BTU/lb/°F) 
Tensile  strength  = 3.09  GPa  (448  ksi) 

Compressive  strength  = 3.09  GPa  (448  ksi) 


TABLE  II:  HMHS  Matrix  Properties: 


Matrix  density  = 3.40E-7  Kg/m3  (0.0457  lb/in3) 
Normal  modulus  = 4.14  GPa  (600  ksi) 

Poisson’s  ratio  =0.34 

Coefficient  of  thermal  expansion  = 0.72E-4/°C 
(0.40E-4  /°F) 

Heat  conductivity  = 0.648  kJ-m/hr/m2/°C 

(0.868E-2  BTU-in /hr /in2 /° F ) 

Heat  capacity  = 1.047  KJ/Kg/°C  (0.25 

BTU/lb/°F) 

Tensile  strength  = 71.0  MPa  (10.3  ksi) 
Compressive  strength  = 423  MPa  (61.3  ksi) 

Shear  strength  = 161  MPa  (23.4  ksi) 

Allowable  tensile  strain  = 0.02 

Allowable  compressive  strain  = 0.05 

Allowable  shear  strain  = 0.04 

Allowable  torsional  strain  = 0.04 

Void  conductivity  = 16.8  J-m/hr/m2/°C  (0.225 

BTU-in/hr/in2/°F) 

Glass  transition  temperature  = 216°C  (420°F) 
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LIST  OF  FIGURE  CAPTIONS  Figure  14  Damage  Progression  under  External  Pressure 

Graphite-epoxy:  AS-4/HMHS[0/90/±4% 

Figure  1 Computational  Simulation  Cycle  Solid  line:  unstiffened  shell 

n ,,  -m*-  ii  *ii  ™ , Dashed  line:  integrally  stiffened  shell 

Figure  2 Finite  Element  Model  with  Square  Elements 

Diagonal  Lines  are  along  the  Integrated  Stiff- 
eners; AS-4/HMHS[0/90/±45]s 
Width  = 305  mm;  Length  = 457  mm 

Figure  3 Deflections  of  Panels  due  to  Bending 

Graphite-epoxy:  AS-4/HMHS[0/90/±45]g 

Solid  line:  flat  panel 

Dashed  line:  integrally  stiffened  panel 

Figure  4 Damage  Progression  under  Uniaxial  Tension 
Graphite-epoxy : AS-4/H M HS [0/90 /±45]S 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 

Figure  5 Displacements  under  Uniaxial  Tension 

Graphite-epoxy:  AS~4/HMHS[0/90/±45]g 

Solid  line:  flat  panel 

Dashed  line:  integrally  stiffened  panel 

Figure  6 Damage  Progression  under  Compression 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]g 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 

Figure  7 Displacements  under  Uniaxial  Compression 
Graphite-epoxy:  AS-4/HMHS[0/90/±4% 

Solid  line:  flat  panel 

Dashed  line:  integrally  stiffened  panel 

Figure  8 Damage  Progression  under  In-plane  Shear 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]g 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 

Figure  9 Displacements  under  In-plane  Simple  Shear 
Graphite-epoxy : AS-4/HM H S[0 /90 /±45]g 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 

Figure  10  Finite  Element  Model  of  Cylindrical  Shell 
AS-4/HMHS[0/90/±45]s 

Figure  11  Damage  Progression  under  Axial  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/db45]g 
Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 

Figure  12  Exhusted  Damage  Energy  under  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 
Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 

Figure  13  Damage  Energy  Release  Rates  under  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 
Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 


Figure  1 Computational  Simulation  Cycle 


Figure  2 Finite  Element  Model  with  Square  Elements 

Diagonal  Lines  are  along  the  Integrated  Stiffeners 
AS-4/HMHS[0/90/±4%;  Width  = 305  mm;  Length  = 457  mm 
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Figure  3 Deflections  of  Panels  due  to  Bending 

Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 

Solid  line:  flat  panel 

Dashed  line:  integrally  stiffened  panel 


TENSION  (kN) 


Figure  4 Damage  Progression  under  Uniaxial  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 


Figure  5 Displacements  under  Uniaxial  Tension 

Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 

Solid  line:  flat  panel 

Dashed  line:  integrally  stiffened  panel 


COMPRESSION  (kN) 


Figure  6 Damage  Progression  under  Compression 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 


DISPLACEMENT  (mm) 


DISPLACEMENT  (in) 
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COMPRESSION  (kN) 


Figure  7 Displacements  under  Uniaxial  Compression 
Graphite-epoxy:  AS-4/HMHS[0/90/±4% 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 


SHEAR  (kN) 


Figure  8 Damage  Progression  under  In-plane  Shear 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 


SHEAR  (kN) 


Figure  9 Displacements  under  In-plane  Simple  Shear 
Graphite-epoxy:  AS-4/HMHS[0/90/d=45]$ 
Solid  line:  flat  panel 
Dashed  line:  integrally  stiffened  panel 


Figure  10  Finite  Element  Model  of  Cylindrical  Shell 
AS-4/HMHS[0/90/db45]s 


DISPLACEMENT  (mm) 
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AXIAL  TENSION  (kN) 
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Figure  11  Damage  Progression  under  Axial  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]s 
Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 


Figure  13  Damage  Energy  Release  Rates  under  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/ ±45]$ 
Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 


AXIAL  TENSION  (kN) 


Figure  12  Exhusted  Damage  Energy  under  Tension 
Graphite-epoxy:  AS-4/HMHS[0/90/±45]S 

Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 


EXTERNAL  PRESSURE  (MPa) 


Figure  14  Damage  Progression  under  External  Pressure 
Graphite-epoxy:  AS-4/HMHS[0/90/ ±45]g 
Solid  line:  unstiffened  shell 
Dashed  line:  integrally  stiffened  shell 


DERR  (kPa) 


